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Abstract 

The generalized dynamics describing the propagation of neutral kaons 
in randomly fluctuating media is derived and analyzed. It takes into ac- 
count possible matter-induced effects leading to loss of phase coherence 
and dissipation. The study of selected neutral kaon observables indicates 
that these non-standard effects are amenable to a direct experimental 
analysis. 
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1. INTRODUCTION 



When a quantum system is immersed in an external environment, its time evolution 
is no longer unitary, since in general it exchanges energy and entropy with the larger 
system. In principle, the subdynamics can be readily obtained from the total, unitary 
system+environment evolution by eliminating {i.e. integrating over) the degrees of freedom 
pertaining to the environment. In practice however, this procedure leads to tractable 
dynamical equations only when the interaction between subsystem and environment can 
be considered weak: in these cases memory effects disappear, and the time evolution of 
the subsystem is described by linear maps, encoding non-hamiltonian physical effects, like 
irreversibility and dissipation. [1-4] 

The system plus environment paradigm for the treatment of the so-called open sys- 
tems is nevertheless very general and has been successfully adopted to model very different 
physical situations, in laser, atomic and molecular physics. [1-8] In particular, it has been 
very useful in describing the effects of random media or of stochastic external fields in 
particle propagation inside interferometric devices. [9-13] Indeed, for weakly-coupled sys- 
tems, the decoherence and dissipative phenomena induced by the media are in general very 
small, so that the most suited way to study them is through appropriate interferometric 
set-ups. [14-16] 

Motivated by these results, in the following similar techniques will be adopted to 
study effects of irreversibility and dissipation in neutral kaons propagating inside a random 
medium. Because of strangeness oscillations, neutral kaons have been one of the prime 
laboratory for the discovery and analysis of small physical effects, CP violation being the 
most striking example; [17,18] therefore, they appear to be the natural place for studying 
environment-induced dissipative effects, also in view of the unique opportunities offered 
by the production of correlated kaons at (/)-factories.[19] 

Matter effects in the physics of neutral kaons have been studied since the very early 
days of kaon physics. However, most of the attention has been devoted to the so called 
regeneration phenomena, where short-lived neutral kaon states are regenerated in a beam 
of purely long-lived ones via the passage through a thin slab of material: this is due 
to the coherent interactions of the kaons with the nuclei of the material, that produce 
different scattering amplitudes for the kaon and antikaon components of the impinging 
beam. One can show that this effect is dominated by forward scattering processes, and 
it can be described by an effective hamiltonian involving the index of refraction of the 
material. [20-23] 

Here, we shall generalize this physical situation by considering the propagation of 
neutral kaons in a randomly fluctuating medium. The kaon system can then be viewed as 
an open system, where the environment (the medium) is described by a classical, random 
external field. 

Quite in general, any environment can be modelled in this way, provided the charac- 
teristic decay time of the associated correlations is sufficiently small with respect to the 
typical evolution time ts of the subsystem. In the case of the kaons, ts can be roughly 
identified with the lifetime of the short-lived kaon, so that the correlations in the material 
through which the kaons propagate must decay very rapidly, or equivalently, the medium 
must fiuctuate on times much shorter than ts- Although this condition looks quite restric- 
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tive, it can be met quite easily by a careful choice of material. Indeed, many short-time 
physical phenomena, like molecular vibrational motion and relaxation, or collisions in liq- 
uids, take place at times that are at least a couple of orders of magnitude smaller than r^.^-' 
Therefore, a kaon moving e.g. in a gas at sufficiently high temperature or in a liquid would 
see a random fluctuating distribution of centers of scattering, and therefore be subjected 
to stochastic incoherent interactions with the material; as a consequence, its dynamics can 
be effectively modelled as being that of a particle propagating in a random medium. 

This physical situation is clearly different from the one encountered in the regeneration 
phenomena, where a single coherent scattering in the thin slab is enough to account for the 
effect. In the present case instead, correlations in the material play a fundamental role: as 
we shall see, they are responsible for the generation of irreversibility and loss of quantum 
coherence. 

In the next section, we shall discuss in detail the derivation of the master equation that 
describes within the weak-coupling hypothesis the propagation of neutral kaons in random 
media. Since it incorporates the presence of dissipative phenomena, it can not be fully 
written in the familiar hamiltonian form; rather, it assumes the structure appropriate for 
generating a quantum dynamical semigroup. [1-4] Specific examples are presented in Sect. 3, 
while some general properties of its solutions are discussed in Sect. 4. Since the matter 
induced effects are small, a suitable perturbative approximation can be used, allowing 
the explicit evaluation of relevant kaon observables: Sect. 5 and 6 will be devoted to the 
analysis of the behaviour of these observables. We shall first study the decays of single 
kaons, relevant for fixed-target experiments, and then discuss the case of correlated kaons 
at (/)-factories. In both cases, the matter- induced phenomena modify in a very specific 
way the various kaon observables, that are therefore amenable to a direct experimental 
analysis. In particular as discussed in the final Sect. 7, at (/)-factories one can adjust the 
experimental conditions so as to allow the propagation in the medium of only one of the 
two correlated kaons; in this configuration, one can easily perform tests on the physical 
consistency of the adopted treatment, shedding further light on the general description of 
open quantum systems. 

Although preliminary, we hope that the results presented in our investigation will 
stimulate further studies on the dynamics of neutral kaons in random media, encouraging 
in particular a detailed experimental analysis. 



^) The direct study of these very short-lived phenomena have been made possible by 
the recent introduction of lasers that are able to produce pulses reaching the femtosecond 
scale; for a review, see [24] . 

Using a more phenomenological approach, non-standard effects in the propagation of 
neutral kaons in a stationary, low-density material have been recently discussed in [25] . 
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2. MASTER EQUATION 

We shall work within the familiar effective description of the neutral kaon system which 
requires the introduction of a two-dimensional Hilbert space; [17-19] the set 
constitutes a convenient basis in this space. With respect to this basis, single kaon states 
can then be represented by density matrices, i?, i.e. by hermitian 2x2 matrices with 
non-negative eigenvalues. Their time evolution is described in terms of an equation of a 
standard Liouville-von Neumann form: 

^ = R{t) + iRit) iy(°)t + [^(^)] . (2.1) 

The first two pieces on the r.h.s. give the usual hamiltonian contribution, while the addi- 
tional linear map L^, explicitly given in (2.10) below, takes into account the presence of 
the stochastic medium. 

The effective hamiltonian H^'^\ the so-called Weisskopf-Wigner hamiltonian, describes 
the propagation of the kaons in vacuum;^) it contains a non-hermitian term, characterized 
by the natural width of the physical states. The entries of if*^*^^ can be expressed in terms 
of its complex eigenvalues, A^"* = m^'' — i^g'^ /2, A^'' = m^'' — ij^^ /2, and the complex 
parameters p^g \ q^^\ p^^\ q^^\ appearing in the corresponding (right) eigenstates, 

1 , 

(2.2) 

1 . 

For later convenience, it is useful to introduce the following positive combinations, involving 
the eigenstates masses and widths, 

Am(«) = mf - mf , Ar(o) = 7^°^ - , F^") = + , (2.3) 

as well as the complex quantities: 

r(o) ^ r(0) ^ .^^(0) ^ = Ar(°) ± 2zAm(°) . (2.4) 

The effective Hamiltonian H^'^'^ can be diagonalized using the similarity transformation 
induced by (2.2): 

^(0) ^^(0)^(0)^(0)-! ^ (2.5) 
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y(0) 



^(0) _.(o) 



^(0) 



Al?) 



(0) 



A^ 



(2.6) 



^) Throughout the paper, we shall append a superscript to all quantities referring to 
neutral kaons in vacuum, i.e. in absence of the medium. 
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Then, one can write: 



(0) 



the two complex parameters 
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(0) 



A 
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.(0) ^(0) 
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SO) , ^(0) 
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r(0) 



^(0) „(0) 



(2.7) 



(2.8) 



involving the ratios, 



(0) 

JO) _ Psl 



rh = 



SO) 



(0) 

Pl 

(0) 



(0)1 



l)/(|(j(o)| + 1) 7^ 0) violating effects 



(2.9) 



signal CPT (for ^(°) 7^ 0) and T (for ^(°) = (|c7 
in mixing, respectively. 

All this is valid in vacuum; when the neutral kaons propagate in matter, the interaction 
with the nuclei of the material gives additional contributions to the evolution equation, the 
term Lt in (2.1). As explained in the Introduction, we shall consider the case of rapidly 
fluctuating media, that can be represented by classical stochastic fields. The action of the 
media on the kaons can then be expressed in generalized hamiltonian form: 



Lt[R{t)\ = -iF{t)R{t)+iR{t) [F{t)\ 



(2.10) 



where 



/i=0 



(2.11) 



while (Jo is the 2x2 unit matrix and cxi, a"2, the Pauli matrices. The quantities F^{t), 
^ = 0, 1, 2, 3, generate a complex, Gaussian stochastic field; they are assumed to have in 
general non-zero mean, but translationally invariant correlations (a star means complex 
conjugation): 



G^^it -s) = (F^it) F,{s)) - (F^it)) {F,{s)) , 
W^.{t -s) = {F^it) F*(s)) - {F^it)) (F*(s)) . 



(2.12a) 
(2.126) 



Since the generalized hamiltonian F{t) in (2.11) involves stochastic variables, the den- 
sity matrix R{t), solution of the total equation of motion (2.1), is also stochastic. Instead, 
we are interested in the behaviour of the reduced density matrix p{t) = {R{t)) which is 
obtained by averaging over the noise; it is p{t) that describes the effective evolution of the 
kaon states in the medium and allows computing the behaviour of relevant observables. 
We shall now explicitly describe the derivation of an effective master equation for p{t), 
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making the additional assumption that kaons and noise be decouple at t = 0, so that the 
initial state is p{0) = {R{0)) = ^(0)>) 

Since the the hamiltonian H^^^ is statistically independent from the stochastic vari- 
ables, one can choose to average over the noise in the interaction representation, where we 
set* 

R{t) = e'*'''°' %)e-^*^^°'' , (2.13) 



so that: 



— ^3 3 

^ = Lt [m] ^'^(^) ^m(^) m + z m ^ F;{t) [a,{t)] ^ , (2.14) 



dt 

IJ,=0 ij.=0 



with cr^(t) — e** ^-ttH rpj^^ tune evolution of the reduced density matrix in 

the interaction representation, p{t) = (i?(t)), can then be expressed as a series expansion 
involving multiple correlations of the operator L^: 



m = M,[p(0)] = J^Mkim] , (2.15) 

fe=0 

Mk[p]= dsi ds2--- dsk {Ls.Ls, ■ ■ ■ Ls^)[p] . (2.16) 
Jo Jo Jo 



The sum Mt in (2.15) can be formally inverted, and a suitable resummation gives (a dot 
represents time derivative): [26] 



dm 

dt 



Mt M,-'[p(t)] = {Mi + (M2 - Ml Ml) + . . . }[p(t)] . (2.17) 



As mentioned before and discussed in more detail at the end of the section, the action 
of the medium on the travelling kaons is weak. Therefore, one can focus the attention 
on the dominant terms of the previous expansions, neglecting all contributions higher 
than the second-order ones. Further, since the characteristic decay time of correlations 
in the medium is by assumption much smaller than the typical time scale of the system, 
the memory effects implicit in (2.17) should not be physically relevant and the use of 
the Markovian approximation justified. This is implemented in practice by extending to 
infinity the upper limit of the integrals appearing in M2 and Mi (compare with (2. 16)). [1-4] 

By returning to the Schrodinger representation, one finally obtains the master equation 
generating the time evolution of the reduced density matrix p(t) = {R{t)). It takes the 
following explicit form: 

^ = -iHp{t)+HAt)H^' + mt)] , (2.18a) 



This condition is always satisfied in a typical experimental situation, where the kaons 
enter the medium after being produced. 
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where 



1 3 ^ r 



(2.186) 
(2.18c) 



The effective hamiltonian in matter H differs from the one in vacuum i?*^^) by first order 
terms (coming from the piece Mi in (2.17)) involving the noise mean values: 



At=0 



(2.19) 



and by second-order contributions (coming from the second-order terms in (2.17)) involving 
the noise correlations (2.12) through the time- independent combinations: 



3 poo 



A=0 
3 /-oo 



dt 



W^x{t) UU-t) + W:x{t) Ux^{-t) 



(2.20a) 
(2.206) 



where the 4x4 matrix U{t) is defined by the following transformation rule: 
Explicitly, one finds: 



(2.21) 



3 |- ^3 

+ ^ iCio -i(Bio + Boi^ + 2 ^ (^'^^'^ ~ iTmCjk^ 



j,k=i 



(2.22) 



On the other hand, the additional piece L[p] in (2.18c) is a time-independent, trace- 
preserving linear map involving the hermitian 3x3 submatrix Cij, obtained from the 
coefficient matrix (2.206) by letting /i, = 1,2,3. It introduces irreversibility, inducing in 
general dissipation and loss of quantum coherence. Altogether, equation (2.18) generates 
a semigroup of linear maps, Ft : /5(0) i— > p{t) = /t[/3(0)], for which composition is defined 
only forward in time: Ft o Fs — Ff+s, with t,s > 0; it is usually referred to as a quantum 
dynamical semigroup. [1-4] 

As a further remark, notice that when correlations in the medium are negligible, 
i.e. the combination in (2.12) are vanishingly small, the physical situation described by 
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equation (2.18) corresponds to that of the regeneration phenomenon. Indeed, in this case 
the presence of matter is signaled solely by the shift H^^^ in the effective hamiltonian. 
For a kaon-medium interaction dominated by coherent forward scattering, H^^^ becomes 
diagonal and its components turn out to be expressible in terms of the forward scattering 
amplitudes fx, fj^ ^^e medium: [20-23] 



^(1) ^ 

rriK 



Ik 
Ik 



(2.23) 



where rriK is the mass, while v represents the density of scattering centers in the 
medium. 

The amplitudes fx, f-^ have been directly measured in experiments for many types of 

materials, [27, 28] and quite consistently one finds: — 10~^ Ar^^-*, where Ar*^°^ ~ 7^*^-* 

is the typical inverse timescale describing the kaon evolution in vacuum, as generated by 
the hamiltonian H^^\ In view of the fact that H^^^^ coincides with the mean value of 
the stochastic hamiltonian in (2.11), one concludes that also the averages (F^) should 
be in modulus of the same order of magnitude. Noise correlations are however much 
smaller: from the definition (2.12), one can safely deduce the following rough estimate: 
|G| ~ \W\ ~ \{F,j,)\'^ ~ 10~^ (Ar^^))^. Therefore, the weak coupling hypothesis adopted 
in deriving the master equation (2.18) appears physically justified. 



3. EXAMPLES 



Before analyzing in detail the evolution in time of the density matrix p{t), we shall 
discuss some physically interesting instances of the master equation (2.18); they correspond 
to specific realizations of the medium through which the kaons propagate. 

Some general considerations apply to all cases. As discussed at the end of the previous 
section, the noise contributions to (2.18) are expected to be small, in particular those 
involving noise correlations. Therefore, whenever the correlations in (2.12) are multiplied 
by other small parameters, e.g. those coming from the hamiltonian H^'^^ in (2.7), one 
can safely neglect them in comparison with the dominant pieces. These considerations 
are particularly relevant in the evaluation of 13^,^ and C^jy in (2.20). These quantities are 
linear in the noise correlations: the entries of the matrix U{t) appearing in (2.20) can 
then be computed at lowest order. Using the diagonalization (2.5) to write e'^*^'^°^ = 

y{0) ^itH^°^ y(0)-l f^^^^j^g^ g^^^^^^g ^ ^ ^(0) ^ ^(0) = (2.6), 5) from 

(2.21) one explicitly finds: 



U{t) = 



(Jo 
a{t) 



a{t) = 



cos ujt sin cot 
— sin cot cos u>t 



(3.1) 



These conditions, that come from the assumption of CPT and T (hence CP) conser- 
vation in mixing, imply a specific phase choice in the definition of the basis state vectors 
in the Hilbert space; in practice, this poses no problems, since physical observables, being 
the result of a trace operation (see later), are manifestly phase-invariant. 
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where uj = iAT^_' /2. 

Further, the stochastic medium fluctuates on time intervals much shorter than the 
typical kaon evolution timescale 1/ AT^'^K Correspondingly, the noise correlations in (2.12) 
can be taken to have an exponentially decreasing form, with decay parameter A much 
larger than AF^^^. This allows neglecting all higher-order contributions in AF^°)/A while 
evaluating explicitly the coefficients B/^i, and C^i^. 



3.1 Generalized regeneration 

In this case, the stochastic hamiltonian in (2.11) is taken to be diagonal, so that only 
the components with = 0, 3 of the stochastic variables -F'ju(^) are non vanishing: 



F+{t) 
F_(t) 



F^it) = ^{Fo{t) ± Fs{t)) . 



(3.2) 



Further, we assume the correlation functions to have the form (/x, = 0, 3): 



Gij.v{t ~ — G^u e 



-Ai|t-s| 



(3.3a) 

(3.36) 



where G^y and W^j^ are time-independent complex matrices, with W = . Physically 
speaking, this case corresponds to a generalization of the familiar regeneration situation. 
The kaons are still mainly forward scattered by the medium; however, while travelling 
in it, they encounter density fluctuations, whose correlations have the general behaviour 
(3.3). 

With the help of (3.1) and (3.3), one can now explicitly compute the noise contribu- 
tions in (2.18). To lowest order in the small parameters, one flnds that only the i = j = 3 
entry of the matrix dj is non vanishing. 



533 = ?H^>0, 



(3.4) 



while the hamiltonian contributions H^^^ and H^'^^ turns out to be diagonal: 



, , Woo + Wss . Goo + G3S 
{Fo) + I ^ I- 



2ttu 
rriK 



fK 





A2 


fK. 



Ai 



(To + 



^"^'^ + XT 



03 



(T3 



(3.5) 

The complex parameters fx and fj^ now contain both the standard forward scattering 
amplitudes, described by (Fq) and (F3), and terms originating from the correlations (3.3). 

The total effective hamiltonian H can be expressed as in (2.7) in terms of its eigen- 
values A5, Xl and the quantities rs, tl, $, a in matter, deflned through its eigenvectors as 
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in (2.2), (2.8) and (2.9). The differences between these parameters and the corresponding 
ones in vacuum can be expressed as a power series expansion in the combination: 



V 



rriK 



Ik - f-K 
(0) _ 



A 



(3.6) 



for many materials, one finds that \r]\is of order IQ-^H-lQ-^.p?, 28] Although a = ct^°) and 
therefore T violation is unaffected by the medium, to lowest order one obtains: 9 = d^^^ — rj, 
so that, not surprisingly, CPT invariance is effectively broken by matter effects even if it 
is preserved in vacuum; in addition, one has: rs = — Vj "i^l = t^l^ + V- 



3.2 Diagonal correlations 

The noise hamiltonian F{t) in (2.11) is no longer diagonal; in this case the kaons un- 
dergo incoherent scatterings not exclusively in the forward direction. However, we assume 
that the only non vanishing correlation functions in the medium be diagonal, and precisely: 



Wii{t -s) = Wooit -s) = Wii e-^^l*-"l , 
W22{t - s) = Wssit - s) = W22 e-^2l*-"l , 
G^^it - s) = G^^ e-^3l*-^l , = 0, 1, 2, 3 . 

The noise contribution to the effective hamiltonian takes now the form: 



(3.7) 



^ A3 



/u=0 



(3.8) 



i=l 



while the matrix dj which characterizes the dissipative term in (2.18) becomes diagonal: 

(3.9) 



Cij — 



2W22 
A2 



2W22 
A2 



3.3 General case 

In general however, the fluctuations in the medium are randomly directed so that the 
correlation functions in (2.12) assume the generic form (3.3), with all components of the 
constant matrices G^i, and W/^i, non vanishing. In this case, keeping again only dominant 
contributions, one finds: 

B,. = ^, C,.= ^-^. (3.10) 

Ai A2 

As clear from the definitions (2.126) and (3.3), the matrix W turns out to be hermitian and 
also positive. The same properties clearly hold for the 3x3 coefficient submatrix Cij char- 
acterizing the dissipative contribution in (2.18) as given in (3.10), as well as for the previ- 
ously discussed cases (3.4) and (3.9). This makes the semigroup Ft : p(0) 1— > p{t) = /^^[^(O)] 
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generated by the master equation (2.18) completely positive. This property is crucial in as- 
suring the consistency of the generalized dynamics Ft in all possible situations; [1-3, 15, 29] 
as we shall see in discussing correlated neutral kaons, lacking of it may lead to physically 
unacceptable results. It is therefore reassuring that in our treatment the property of com- 
pletely positivity emerges naturally, without further assumptions. ^-^ 

Although the matrix Cij turns out to be in general complex, in the following we shall 
limit our considerations only to situations for which Cij results real. These correspond to 
the physically most interesting cases; in fact, the reality condition guarantees the increase 
of the von Neumann entropy, dS/dt > 0, 5 = — Tr(pln p),[29] which is a physically desirable 
requirement for the dynamics of a small system, the neutral kaon, in weak interaction with 
a large environment, the medium. For later convenience, we shall parametrize the entries 
of Cij using the real constants 



(3.11) 



a 


= Cn 


+ C22 


h = 


C23 


a 


= Cii 


+ C33 


c = 


-Ci3 


7 


= C22 


+ C33 


P = 


C12 . 



The positivity condition on Cij then requires: 



2R = a + j-a>0, RS-b'^ >0 , 

2S = a + j - a>0 , RT - > , 

2T = a + a- 7>0, ST - P'^ > , 

RST -2bcP- RI3^ - Sc^ -Tb^>Q . 



(3.12) 



3.4 White noise 

Finally, let us consider the case of a medium that behaves as a white noise; it is 
described by stochastic variables Ff^{t) that are 5-correlated in time: 

G^u{t - s) ^ G^u 5{t - s) , (3.13a) 
W^,{t-s) = W^^. 5{t-s) , (3.136) 

with G^i/, W^jy time-independent matrices, such that W = . From the definition (2.20), 
one immediately finds: 

B^u = G^, , C^, = 2W^, , (3.14) 

However, when the noise correlations assume the more general exponential behaviour 
Wij,^{t — s) — Wfj^i, e"^'^"'*"^', Apjy > 0, the corresponding matrix Cij is no longer auto- 
matically positive. In this case, complete positivity needs to be imposed as an additional 
requirement. 
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without any approximation, and all considerations about the properties of the generalized 

dynamics Ft presented in the previous subsection can be repeated here. The present case 
is nevertheless special: because of the correlations (3.13), all higher order terms in the 
expansion (2.17) identically vanish, so that the evolution equation (2.18) is in this case 
exact. 



4. OBSERVABLES 

The evolution in time of the neutral kaons in the medium is fully described by the 
reduced density matrix pit), solution of the master equation (2.18): any physical property 
of the system can be extracted from it by taking its trace with suitable hermitian operators. 
Of particular interest are the observables Of that are associated with the decay of neutral 
kaons into a final state /, typically 27r, Btt and semileptonic ttIv states. In the fixed 
\K^)} basis. Of can be represented by a 2 x 2 matrix, whose entries are expressed 
in terms of the two independent decay amplitudes in vacuum A{K^ /), A{K^ — > /):^^ 



Of = 



|^(KO^/)P [A{K^ ^ f)X A{K^ ^ f) 

A{K^ ^ /) [A(K^ ^ /)] * \A(K^ ^ /)p 



(4.1) 



The corresponding physical observable, directly accessible to the experiment, is given by 
its mean value: 

{df)it) = Tr[dfp{t)] , (4.2) 

whose evolution in time is regulated by that of p{t). 

Explicit solutions p{t) of the master equation (2.18) can be obtained using perturba- 
tion theory. As discussed in the previous section, the entries of the coefficient matrix Cij 
characterizing the dissipative contribution L turn out to be of the form Wij/X, where Wij 
are the coefficients of the noise correlations, while A represents their typical inverse decay 
time. From the analysis at the end of Section 2, it follows that their magnitude can be 
evaluated to be of order: \Cij \ ~ 10"'* Ar*^°^(Ar*^'^)/A), which is small, since by hypothesis 
A ^ Ar*^°). These considerations allow treating the dissipative piece L in (2.18) as a per- 
turbation to the contributions of the effective hamiltonian H, which contains the standard 
Weisskopf-Wigner term H^^^ of (2.7), besides the smaller, higher order noise contributions 
ii-W and 

It is convenient to work in a basis in which H is diagonal: all physical quantities (4.2), 
being the result of a trace operation, are basis independent. This can be obtained by a 
similarity transformation analog to the one in (2.5): 

H = VHdV-^, (4.3) 



In the adopted treatment, kaon decays, being the result of weak interaction processes, 
are unaffected by the presence of the media; all media influences are encoded in the dy- 
namical equation (2.18). 
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with 



V = 



PS 
QS 



PL 
-QL 











(4.4) 



where, as before, Xs — ms — ^75/2, Xl = m-L — ^7l/2 are the two eigenvalues of if, while 
PS, PL, Qs, Ql define the corresponding eigenvectors \Ks), \Kl) through definitions similar 
to those in (2.2). 

This change of basis induces the transformation p ^ p = V~^pV~^^ on the density- 
matrix, while the master equation (2.18) becomes: 



dpjt) 
dt 



-iHDp{t) + ip{t)Hj, + L[p{t)] , 



(4.5) 



where L[p] is again of the form (2.18c), but with a new coefficient matrix C^, whose entries 
are linear combinations of the ones of the original matrix Cij {cf. (3.11)).^^ 

Similarly, also the observables Of get transformed. Of ^Of — V'^OfV, so that the 
corresponding mean value (Of) remains indeed unchanged: 



(Of) ^ (Of) = Tt Pit) Of 



(4.6) 



The entries of the transformed matrix Of can be conveniently expressed in terms of the 
following two complex quantities: [17, 18] 



f _ qs A{K^ ^ /) 
PS A{K^ ^ /) ' 



as 



Of^\A{K'^f)\' 



\ps?\l + X 



/|2 



_ QL AjRO f) 
^ PL A{K^ ^ /) ' 



PsPL (1 + A^)*(1-A{) 



(4.7) 



.PsPl{l + Xi){l-X{y \PL?\l-Xi\^ 



(4.8) 



or alternatively in terms of their reciprocal fig = ^/Xg, = 1/A{^, so that 
Of = \A(K^^f)f 



qs\^\l+pi\^ -qs<lL{l + piy {I-pD 



-qsql{l + pi){l-piy \qL\^\l-fii\ 



(4.9) 



Using standard perturbation theory, solutions of (4.5) can then be easily obtained to 
any needed order in the small parameters Cij, while keeping an exact dependence on the 
quantities parametrizing the effective hamiltonian H, i.e. masses ms, ttll and widths 75, 



Explicitly, one has: L[-] = y L[y ■ 1/^] V"^; however, since the coefficients Cij in L 
are small, in keeping with the approximation used before, one can neglect CP and CPT 
violating effects in matter, and set ps = Pl = qs = Ql = l/v^ in the above transformation 
law. 
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7l in matter, and the ratios rs = Ps/qSi i^l = Pl/ql- While these constants describe 
famihar hamiltonian behaviour plus CPT and T violating effects in matter [through the 
combinations 6 and ^ = {\a\ — l)/{\a\ + 1)], the additional piece L in (4.5) is responsible 
for matter-induced non-standard effects, leading to dissipation and quantum decoherence. 
In the following we shall focus on the latter, analyzing in detail how the new effects modify 
the behaviour of selected neutral kaon observables. 



5. SINGLE NEUTRAL KAONS 



We shall start by discussing the time evolution of observables associated with the 
propagation in random matter of single, uncorrelated neutral kaons; typically, these can 
be studied in fixed target experiments {e.g. see [30]). As discussed in the previous section, 
relevant physical quantities are the probability rates that an initial neutral kaon decays at 
a certain time t into a given final state /: 



(5.1) 



where pK°{t)i PKo{t) represent the evolution according to (4.5) of initially pure 
mesons states. 

The case of semileptonic final states is of particular interest. The amplitudes for the 
decay of a or state into tv~£~^u, Tv~^£~iy can be parametrized as follows: [31] 



A(K^ 7T+e-i?) ^M%1 + y*) , 



(5.2a) 
(5.26) 
(5.2c) 
(5.2(i) 



where At is a common factor, while the parameters a;, z measure violations of the AS' = AQ 
rule and y signals violations of the CPT-symmetry. These quantities are expected to be 

very small, so that one can consistently neglect terms containing x, y ox z when multiplied 
by other small parameters, coming either from the hamiltonian {e.g. 9 and ^) or the 
dissipative part of the evolution equation (4.5). In particular, this approximation implies 



that for semileptonic decays: A^^ = 

Then to first order in the dissipative matter effects, the probabilities (5.1) explicitly 
read:^^ 



^1 



\M\' 



e ^^cos(Amt) 



i^^^^e-(^-^)^^*-27ee(y + 2C) 



-he~^* sin(Amt) 



AXm{rs r£) 



2Im{x) +ne{B) 



Since AF ~ 2Am to a good approximation, in writing (5.3) and (5.4) we have set the 
ratio 2Am/Ar equal to one when multiplying small parameters. 
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+e 



+e" 



-7-L* 



-1st 



- ne{x + y-2C) + D 



+ ne{x -y + 2C)-D 



(5.3a) 



\M\' 



e ^*cos{Amt) 



+e ^*sin(Amt) 



+e" 



+e 



-7st 



2 



4g-(A-D)Art 



2Xm(2 + 2C) -7ee(S) 

+ 7^e(y - - 2Tm(C) + D 
2 + 7^e(y + + 2Xm(C) - D 



(5.36) 



where Am, AF and F are defined as in (2.3) in terms of masses and widths in medium. 
Matter induced dissipative effects are controlled by the dimensionless parameters: 



a + a 



B 



a — a + 2ib 
Am 



C 



c + i(3 



D 



7 

AF 



(5.4) 



The expressions for Vi+{K^;t) and V£~{K^;t) can be obtained from (5.36) and (5.3o) 
respectively, by changing the sign of y and C, and letting rs ^/rs, tl ^/fL-, x ^ z. 
For a non fluctuating medium, one has A = B = C = D = and the expressions in (5.3) 
reduce to the standard ones, giving the probability of a semileptonic decay for a kaon that 
has travelled in a slab of material. 

The probabilities Vi are directly accessible to the experiment. Therefore with a suit- 
able set-up, thanks to the different time dependence in the various pieces of (5.3), it is 
possible to extract information on the parameters (5.4) and therefore on matter induced 
decoherence effects. Clearly, this task becomes easier in the case of the generalized regen- 
eration discussed in Sect.3.1: in this case, the stochastic medium fluctuations are such that 
A = D = S/2, C = 0. 

Similar considerations apply to the study of the decay of a single neutral kaon into 2it 
or Stt final states. For instance, in the just mentioned simplified case, the 27r-decay rate is 
sensible to the surviving dissipative constant A: 



''Ty[02^Pko{Q)\ 



\rs 



|2|1 + A|-| 



\rs + rL] 



A 



+ e 



Ws + '^lP 



+ A 



(5.5) 



+ 2e |?727r| cos (Am t - 0277 ) 
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where 



(5.6) 



while A|^, Xj^ are the decay parameters defined in (4.7). 

More in general, combining the semileptonic and pion decays probabilities in suitable 
asymmetries one can obtain enough independent observables to treat the case of more 
general media, for which the parameters in (5.4) are all different from zero. 



6. CORRELATED NEUTRAL KAONS 

Matter induced dissipative effects can be further studied in experiments at (^factories, 
using correlated neutral kaons.[19] Indeed, these set-ups behave like quantum interferom- 
eters and therefore are particularly suitable for analyzing phenomena leading to loss of 
quantum coherence. 

In a 0-factory, correlated kaons are produced from the decay of the meson. Since 
this is a spin-1 particle, its decay into two spinless bosons produces an antisymmetric 
spatial state. In the (j) rest frame, the two neutral kaons are then seen flying apart with 
opposite momenta, and in the basis \K^), \K^), the resulting state can be described by: 

\iPa) = ® \K^,p) - \K^,-p) ® li^^p)) . (6.1) 

The corresponding density operator pA — li^A) (V'^l is represented by a 4 x 4 matrix, since 
now it describes two kaons. The time evolution of the correlated two kaon system can 
be expressed in terms of the single meson dynamics Ft generated by the equation (2.18). 
Indeed, once produced in a 0-decay, the two kaons evolve independently, so that the density 
matrix that describes the situation in which the first kaon has evolved up to (proper) time 
ti, while the second one up to (proper) time t2 is given by: 

PA{tut2) = (^rt,(g)rt^)[pA] . (6.2) 

Correspondingly, one can now study double decay observables, i.e. the probability that a 
kaon decays into the final state /i at time ti, while the other kaon decays at time t2 into 
the final state /2: 

V{h,h;f2,t2) = Tr[(df^0df,) PA{ti,t2)\ ; (6.3) 

here, O/^, Of^ are the 2x2 hermitian matrices introduced in (4.1) that describe the decay 
of a single kaon into the final states /i, /2, respectively. 



10) -pQj- ^]^g actual computation of the probabilities V{fi,ti; /2, ^2), it is again convenient 
to work in the basis in which the effective hamiltonian H is diagonal, and therefore use 
PA = [V-^ ® F-i] PA [W-i yt-i], and Of in (4.8), (4.9). 
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The probability rates in (6.3) are very sensitive to matter induced decoherence effects. 
This is most strikingly shown by considering correlated decays at equal time ti = t2 = t 
into the same final state /i = /2 = /• In absence of the dissipative term L in the evolution 
equation (2.18), the antisymmetry properties of the initial state pA would be preserved by 
the factorized evolution Ft ® Ft, thus producing a vanishing result for V{f ,t; f ,t). The 
equal time probabilities P(/, /, t) are therefore particularly suited to signal the presence 
of the dissipative parameters in (5.4). For instance, for / either a 27r or a Stt final state, with 
associate intrinsic CP parity C/, the probability V{f,t;f,t) is sensitive to the dissipative 
parameter D; explicitly, one finds: 

V{f, t; /, t) oc 2D (e-^-* - e'^^') [(1 + Cf)e-^'' + (1 - Cf)e-^''] , (6-4) 

the proportionality constant being a decay amplitude normalization factor dependenig on 
whether / = 27r or Stt. On the other hand, in the case of semileptonic final states, it is 
the parameter a that determines the slope at which the corresponding joint probability 
approaches zero for small times; indeed, one has: 

V{e'^,t;e^,t) at . (6.5) 

Although these probabilities together with the more general ones in (6.3) can be 
measured at a (/)-factory, much of the experimental analysis performed at these set-ups is 
devoted to the study of integrated distributions at fixed interval t = ti — t2:[S2] 

nOO 

r{h,f2;t)= dt'V{h,t' + t;f2,t') , t>0. (6.6) 
Jo 

A particularly interesting observable that can be constructed with these integrated prob- 
abilities involves 2it final states: 

, _ r(7r+7r-, 27r°; t) - r(27r", tt+tt"; t) 

' r(7r+7r-,27rO;t)+r(27rO,7r+7r-;t) ' ^ ' 

it allows determining the ratio e' /e, where e and e' are the familiar phenomenological 
constants parametrizing the decay in two pions of the short and long-lived kaons in vac- 
uum. [17, 18] In a medium, the asymmetry As>{t) gets new contributions, both from the 
effective hamiltonian H and the dissipative term of (2.18). To lowest order in all small 
parameters, one explicitly finds: 



where 



^'■w--(i)t;i--'"(7)ig. 

N2{t) = 2 |£|^e~^*sin(Amt) , (g g) 

Nsit) = e"^^* (|£p +D)+ e-Ts* (^\e\^ - - 2 \e\^ e'^^ cos(Amt) , 
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and e = e^EL- with sl = {vl - l)/{rL + 1), = in' ' ^)/in' + !)• A 
careful analysis of the time behaviour of the two contributions in (6.8) would provide a 
way to estimate the dissipative parameter D, together with real and imaginary part of 
e' /e. Further, note that in the long time limit the asymmetry (6.8) reduces to: 

/ I ~|2 

A'(r)~37^e(|) , (6-10) 

and not simply to 37le[s' /e) , as in vacuum. Therefore, even assuming El — Sl^, a measure 
of As> can no longer provide a determination of TZe{e' /e) unless an estimate on the matter 
induced dissipative parameter D is independently given. 



7. DISCUSSION 

Neutral kaons propagating in a stochastically fluctuating medium can be treated as 
an open system, i.e. as a subsystem immersed in an external environment. Starting from a 
microscopic hamiltonian with a generic, linear kaon-matter interaction term, a generalized 
subdynamics for the kaon states has been explicitly derived by averaging over the matter 
noise. It takes the form of a completely positive quantum dynamical semigroup, where the 
presence of the medium manifest itself through the generation of hamiltonian corrections, 
that modify the familiar Weisskopf-Wigner description of the neutral kaon system, and 
ii) the addition of extra pieces inducing dissipation and loss of quantum coherence. 

Some of the hamiltonian contributions have been analyzed before in connection with 
the so-called regeneration phenomena: they arise because of the coherent interaction of the 
travelling kaons with the scattering centers in the medium. On the contrary, the remaining 
hamiltonian pieces and the new dissipative contributions in the kaon evolution equation 
originate from the stochastic correlations in the medium: as they move in the material, 
the kaons encounter density fluctuations, whose correlations decay in time very rapidly 
with respect to the typical time scale of the kaon system, thus inducing irreversibility and 
decoherence in their dynamics. 

Many physical phenomena can give rise to short time correlations in ordinary ma- 
terials: they have been studied by the so-called femtosecond chemistry. [24] By suitably 
inserting one of these materials in any standard kaon-physics set-up, one can experimen- 
tally study the new, matter induced dissipative effects. Indeed, as discussed in Sect. 5 
and 6, one finds that both single and correlated kaon decay observables get modified in 
a very specific way by the presence of these effects; as a result, they can be probed quite 
independently from other kaon physics phenomena. 

Although here derived in a specific context, the generalized evolution equation (2.18) 
has wider validity: it has been shown to generate the most general open system dy- 
namics compatible with a semigroup composition law and the requirement of complete 

A non vanishing D would decrease As', making 7le{e' /e) bigger than measured. In the 
case of string induced dissipative effects, this phenomenon have been discussed in detail in 
[33] . Similar conclusions have also been mentioned in [25] . 
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positivity. [1-3] As such, it has been recently apphed to the description of dissipative effects 
induced at low energies by the dynamics of fundamental objects (strings and branes) at a 
very high scale, typically the Planck mass. [34] Although very small, these string induced 
decoherence effects might be experimentally studied using interferometric devices, like 0- 
factories,[35, 36] and can in principle interfere with the phenomena described here. Notice 
however that the two situations correspond to totally different experimental conditions. 
Matter generated phenomena are completely under the experimental control; the effects 
they induce can easily be isolated from those that might be generated at Planck scale by 
suitably varying the experimental conditions. 

The possibility of choosing and tuning at will the experimental set-up further allows 
performing interesting tests on the physical consistency of the dynamics generated by 
(2.18), and in particular on the property of complete positivity. [29] Consider the case 
in which only one of the two correlated kaons coming from a ^-meson decay actually 
propagates in a stochastic medium, while the other evolve in vacuum. The density matrix 
that describes this situation at time t is given by pAit^t) = [Tt ® rf^^)[pA]-, where Ft is 

the map generated by the equation (2.18) evolving the kaon in the medium, while r^'^\ 
generated by the Weisskopf-Wigner hamiltonian H^^^ of (2.7), describes the propagation 
in vacuum of the second kaon. 

As mentioned in Sect. 2, any density matrix must be positive; this requirement comes 
from the physical interpretation of its eigenvalues as probabilities, that thus must be non- 
negative. In the case of a medium with diagonal correlations as discussed in Sect. 3. 2, for 
which the parameters in (3.11) are such that a = a, h = c = j3 = 0, the four eigenvalues 
Ai(t), z = 1, 2, 3, 4, of the matrix p^(t, t) above can be explicitly computed: 



Although Ai(t), A2(t), ^2,{t) are manifestly positive for any t, because of the minus sign in 
front of the square root, one can check that X^it) is non- negative only for 



This is precisely the inequality that in this case guarantees the condition of complete 
positivity of the dissipative evolution Ft (compare with (3.12)); lacking of it would have 
led to physically inconsistent dynamics. 

The situation just described can certainly be realized at a (/)-factory, so that, at least in 
principle, the time behaviour of the above eigenvalues can be experimentally studied, and 
the inequality (7.3) probed. This would allow a direct test of the condition of complete 
positivity, thus providing direct experimental support for one of the crucial properties 
characterizing the quantum dynamics of open systems. 



where 




(7.1) 



(7.2) 



2q; - 7 > . 



(7.3) 
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